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Abstract 

We use temperature tuning to control signal propagation in simple one- 
dimensional arrays of masses connected by hard anharmonic springs and with 
no local potentials. In our numerical model a sustained signal is applied at 
one site of a chain immersed in a thermal environment and the signal-to-noise 
ratio is measured at each oscillator. We show that raising the temperature 
can lead to enhanced signal propagation along the chain, resulting in thermal 
resonance effects akin to the resonance observed in arrays of bistable systems. 
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I. INTRODUCTION 



In the past few years it has become abundantly clear that the presence of noise in non- 
linear systems may lead to an enhancement of a number of often desirable features such as 
energy localization and mobility and the detection and propagation of weak signals. The in- 
terplay of stochasticity and nonlinearity that amplifies the system response is a cooperative 
phenomenon whose detailed nature depends on the particular structure of the system and the 
forces acting upon it [0,0. One manifestation of the interplay is found in the phenomenon 
called stochastic resonance, which has been invoked in a wide range of physical , chem- 
ical |§-|^, geological |I|,pi][, and biological [|,|n|] systems. Recent literature, including our 



own work, [|T2| -[l^ has focused on spatially extended systems [|T^ including noise-enhanced 
propagation in coupled arrays of bistable units 0,0, excitable media ||7|,|5|JT7| , reaction- 
diffusion systems and dynamics and signal propagation in cardiac tissue [0,0 . It has 
been repeatedly noted that discrete extended systems pose particular mathematical chal- 
lenges that have barely been explored in spite of the fact that many physical systems are 
intrinsically discrete |pTl-p5[]. 

The ubiquitous picture of stochastic resonance involves a particle moving in a double- 
well potential subject to a weak external signal that periodically changes the potential by 
alternately raising and lowering the wells The signal is "weak" if the periodic force 
is too small to cause the particle to scale the barrier between the wells. Nevertheless, an 
appropriate random force is sufficient to cause the particle to cross over the barrier even 
in the absence of a deterministic signal. In the simultaneous presence of a weak signal 
and a sufficiently weak noise, the transitions over the barrier occur rarely and at a rate 
determined by the noise intensity. These transitions are slow compared to the frequency 
of the deterministic signal; the transition rate then carries little information about the 
signal. At the other extreme, when the noise is strong it induces rapid transitions that are 
again essentially uninfluenced by the frequency of the signal. At an optimal noise intensity, 
however, the mean first passage time associated with the noise and the frequency of the 
signal are in synchrony (stochastic resonance), and the passage from one well to the other 
carries maximal information about the signal frequency. 

A less ubiquitous but nonetheless important occurrence of stochastic resonance (that 
has been called "nonconventional" by its discoverers) arises for particles moving in nonlinear 
monostable potentials |0j2^ . It is argued that stochastic resonance can be expected to occur 
in any single-well underdamped system for which the spectral density of the fluctuations of 
the system in the absence of a periodic signal exhibits a well-resolved narrow peak that grows 
faster than quadratically with temperature. The effect is conflrmed via analog simulations 
of a single- well Duffing oscillator [ PB||27| and of a SQUID loop More recently, stochastic 
resonance at higher harmonics in monostable systems was ascertained for an overdamped 
system when the nonlinearity is not concentrated at the equilibrium position [^ . 

Recent developments in the field have generalized these ideas to linearly coupled arrays 
of bistable oscillators P,|^JT(3[]. A signal with the help of the noise in these arrays can cause 
a "phase jump." If the noise is sufficiently weak, the phase jump travels in the form of a 
moving kink (strong noise causes random phase jumps that make it difficult to separately 
identify a phase jump associated with the signal). The creation or destruction of a phase 
kink is an activated process, i.e., the signal and/or noise must be sufficiently strong to cause 
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a transition from one well of one of the bistable potentials to the other. The presence of 
such kinks is associated with an "energy gap:" it takes a finite amount of energy to destroy 
a kink. The language used in this description was originally borrowed from the kink soliton 
context. 

Reported instances of stochastic enhancement and stochastic resonance in extended ar- 
rays involve coupled (overdamped) bistable units. Herein we show that enhanced propagation 
can be achieved through thermal tuning of even simpler discrete arrays of masses connected 
by monostable anharmonic springs (with no local potentials). The signal is identified with an 
amplitude that exceeds (by a predetermined amount) that due to the thermal background. 
Here there is no activation process and no energy gap, and the signal can simply disperse or 
dissipate. The language appropriate to this case is akin to that originally associated with 
envelope solitons. We focus on the propagation distance and amplitude along the chain of a 
signal continuously applied at one site of a ID array, and compare results for harmonic and 
hard anharmonic chains. In particular, we show that when the anharmonic chain is immersed 
in a thermal bath, it is possible to maximize the distance of propagation and the amplitude 
of the signal at a given site by tuning the temperature to particular optimal values. We 
call this phenomenon thermal resonance. Our systems are in general not overdamped and 
thus include inertial contributions to the motions of the masses. Noise and damping repre- 
sent a realistic thermal environment with a tunable temperature and dissipation that obey 
an appropriate fluctuation-dissipation relation. We have found and reported elsewhere 



that the propagation of an energy pulse in a hard anharmonic array can be enhanced by 
immersion in a thermal bath, and that hard anharmonicity in the springs causes a tight and 
persistent packing of the energy. Those results suggest the possibility of a thermal resonance 
in the transport of a sustained external signal in these simple arrays. 

In Section |I| we present our model and some details of the numerical integration of the 
equations of motion. Our characterization of a thermal resonance is presented in Section |T|, 
and our main results are shown in Section |rV|. In Section |V| we discuss the dependence of our 
results on different parameter models. Section IVl contains our summary and conclusions. 



II. THE MODEL AND NUMERICAL PROCEDURE 

Our model consists of a one- dimensional chain of unit-mass sites, each connected to 
its nearest neighbors by either harmonic (quadratic) or hard anharmonic (quartic) springs. 
Accordingly, the Hamiltonian of the array is 

^ = E (y + f - ^n-lY + -^{Xn ' a:„_l)^j , (l) 

where k and k' are the harmonic and anharmonic coupling constants respectively. Thermal- 
ization of the chains is achieved through a Langevin prescription for coupling the system 
to a heat bath. The stochastic equations of motion for sites n = 1, ■ ■ ■ , — 1 are obtained 
from the Hamiltonian augmented by the usual Langevin forces: 

Xn k{Xfi—\ Xri) kiyXfi Xn+l) ~\~ k {Xfi—\ S^n) 

- k'{Xn - Xn+lf - IXn + /n(i) , (2) 
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where a dot represents a derivative with respect to time. The fn{t) are zero-centered, Gaus- 
sian, ^-correlated fluctuations that satisfy the fluctuation-dissipation relation at temperature 

{fn{t)fn'{t')) = 2^kBTSn,n'S{t - t') (3) 

{kB is Boltzmann's constant). We impose periodic boundary conditions, so that xn = Xq. 
A sustained signal is applied to the site n = that determines its velocity at all times: 

xq = A sin(co'o^) • (4) 

The positions and momenta of all the other sites are otherwise "free" and determined by the 
equations of motion. We study the propagation of this signal along the chain as a function 
of the temperature. This particular way of applying a signal is of course not unique (e.g. 
one might apply an oscillatory force instead), but we have ascertained that the results are 
insensitive to the detailed choice. 

An analytic solution of this problem is not available for an anharmonic chain, so we must 
rely on numerical integration, which is performed using the second order Heun's method 
(equivalent to a second order Runge-Kutta integration) [^^. The time step is deter- 



mined by the period of oscillation of the velocity of the first site, r = 2tt/ujo, through the 
relation At = For each simulation, the system is initially allowed to relax to thermal 

equilibrium. For all the simulations presented, this is achieved in less than 20 units of the 
dimensionless time. Typically, after a transient that is longer the further the site is from 
site (and thus a measure of the velocity of propagation), each site settles into stationary 
behavior that is a mixture of thermal motion and response to the signal. Sites that are far 
from the signal never exhibit this transient (thus indicating a finite distance of propagation) 
and simply continue their thermal motion. At any given site that is reached by the signal, 
one can observe the amplitude of the motion associated with the signal over and above the 
thermal motion. 

Our interest here lies in demonstrating resonances in the dependence on temperature 
of the propagation distance and velocity and of the response amplitudes once stationarity 
has been achieved (in all our simulations, each site has settled into its long-time behavior 
after 100 first-site oscillations). Possible energy return effects around the periodic chain are 
prevented by making the chain sufficiently long and/or sufficiently increasing the dissipation 
parameter of distant sites. Our chains typically consist of 70 sites with a large dissipation at 
sites 27-32 (these numbers can easily be varied). Our "measurements" are then taken over 
80 oscillation periods and non-zero-temperature results are averaged over 500 realizations. 
All of these choices (equilibration time before applying the signal, integration time step, 
transients, length of chain, and number of realizations) have been carefully tested. 



III. CHARACTERIZATION OF THERMAL RESONANCE 

We must choose a sensible response variable to characterize the behavior of our array. 
When stochastic resonance is studied in arrays of bistable potentials, the system response 
is usually analyzed in terms of a crossover time series that characterizes the transitions of 
each bistable element from one well to the other. In excitable media a reasonable response 
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variable involves firing times of the individual elements. In our system the most convenient 
choice is the velocity of each site because it corresponds directly to the applied signal, and 
because the time average of the velocity in the stationary state vanishes at any temperature. 
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FIGURES 




FIG. 1. Mean velocity (v) of the 12*^ site of a hard anharmonic chain with /c' = 3, 7 = 0.5, 
ujQ = 0.5, and A = 0.5, for different values of the temperature. Enhancement of signal propagation 
is observed with increasing temperature (1** and 2"*^ panels), but the signal is destroyed as the 
temperature further increases (3'''^ and 4*'^ panels). Note the different vertical scales. 

Fig. ID is a dramatic but typical demonstration of a realization of thermal resonance. 
It shows the mean velocity (v) (averaged over realizations) as a function of time for the 
12*'^ site of a hard anharmonic chain at different temperatures. The first panel shows the 
results for zero temperature and for a very low temperature; the temperature increases in 
subsequent panels. At zero temperature the 12^^ site hardly moves because the signal has 
been dissipated to the bath before reaching this site (in the corresponding harmonic chain 
the signal reaches the 12*^ site quite vigorously at zero temperature, a confirmation of the 
fact that a given dissipation is much more effective in a hard anharmonic potential than 
in a harmonic one [lT2| , p!3| , |32| - see Fig. A very small temperature increase (still first 
panel) causes a large enhancement of the signal, which clearly now reaches the site. This is 
apparent in the oscillatory behavior of the velocity over and above the noisy background. 
The temperature in the second panel is close to its optimal value, that is, the value that 
most enhances the signal at this particular site relative to the thermal background. Hence 
the motion of the 12*'' site at this temperature is mostly driven by the periodic forcing of the 
first site. A further increase in the temperature (third panel) causes the average velocity to 



6 



become increasingly noisy because ever larger fluctuations dominate the dynamics. Finally, 
at a sufficiently high temperature (fourth panel), the signal is essentially buried in the 
fluctuations, and the motion is simply that imposed by the thermal bath. In contrast, in a 
harmonic array [k = 3) the signal at any site simply degrades with increasing temperature. 
This is illustrated in Fig. 0. 
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FIG. 2. Mean velocity (v) of the 12^^ site of the harmonic chain with /c = 3, 7 = 0.5, ujq = 0.5, 
and A = 0.5, for different values of the temperature. Signal degradation with increasing tempera- 
ture is clearly observed. Note the different vertical scales. 

To provide a quantitative measure of the thermal resonance, we define the power spectral 
density Sj^uj) at each site j as 



Sj{uj) 



+00 



{vj{t)vj{t + t)) dr 



(5) 



where the brackets denote an ensemble average over realizations and an average over time. 
In addition to the thermal fiuctuations, this function contains the spectral information about 
that part of the signal that has reached site j. An example of a portion of the spectrum 
for the case that we will call our "standard case" {k' = 5, 7 = 0.2, c^o = 1.0, A = 0.5) is 
shown in Fig. |^. The signal extraction from background noise that characterizes stochastic 
resonance is traditionally performed via a signal-to-noise ratio (SNR) p[: 
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(6) 



where the signal power is the value Sj^uq) — Sj^noisei^o), the thermal power Sj^noisei^o) is 
estimated by performing a 4*^-order polynomial fit to Sj{uj) around - but not including - 
the forcing frequency Uq, and Au denotes the frequency integration step and is equal to 
0.0125 (the inverse of the 80 oscillation periods used as our measurement time) throughout 
this paper. This definition of the SNR is not unique, but our results are robust with respect 
to variations in this definition. 
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FIG. 3. Power spectral density at the IS*'^ site of the hard anharmonic chain with k' = 5, 
7 = 0.2, luq = 1.0 and A = 0.5 (we call this our "standard case" in the text) at temperature 
ksT = 0.5. The wide solid line shows the polynomial fitting around o^o- The circle indicates the 
value of 5i5(ti;o) and the square that of Si^^noisei^o)- 

IV. THERMAL RESONANCES 

We first present SNR results for the harmonic chain, so as to clarify later the ways 
in which the anharmonic chain behaves differently. The analytic calculations associated 
with the harmonic chain are presented in Appendix |^. These results serve as a test for 
our numerical simulations. The first panel in Fig. ^ shows SNR curves as a function of 
temperature for different sites, and the second panel shows the same results as a function 
of distance from the forced site for different temperatures. The results are exactly as shown 
in the Appendix and as one would expect: the SNR decreases monotonically with increas- 
ing temperature and with increasing distance from the applied signal. Note that SNR{j) 
decreases with increasing temperature because the numerator in Eq. (^ is essentially inde- 
pendent of temperature while the denominator increases (see Appendix A) . A point to note 
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is that the decay of SNR{j) with j at a given temperature provides a measure of the shape 
of the stationary front of the signal at that temperature. 
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FIG. 4. First panel: Typical SNR curves as a function of temperature for different sites along a 
harmonic chain. The SNR for each site decreases monotonically with temperature. This particular 
example shows sites from the 5*^ to the 14*'' (top to bottom) with fc = 3, 7 = 0.5, ujq = 1.0, and 
A = 0.5. Second panel: SNR as a function of site for different temperatures. 





FIG. 5. First panel: SNR curves for different sites (from j = 6 to 15, top to bottom) along the 
anharmonic chain as a function of the temperature for the standard case: A = 0.5, k' = 5, j = 0.2, 
loq = 1. The optimal temperature increases with distance from the first site. Second panel: SNR 
as a function of the site for different temperatures. 

We now turn to the anharmonic chain. The first panel in Fig. ^ shows SNR curves 
as a function of temperature for different sites. The thermal resonance is identified with 
the j-dependent maximum of SNR{j) as a function of temperature. The signal at the first 
few sites is monotonically weakened as a function of the temperature (as it is at all sites 
in the harmonic chain), but the SNR is enhanced with increasing temperature for sites 
farther away until it reaches a maximum (thermal resonance), beyond which it decreases. 
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Note that the optimal temperature increases with distance from the first site. For the 
particular parameters used in this simulation the thermal resonance occurs most clearly at 
temperatures in the range ksT < 4 for sites between the 8^^ and the 15*''; these details can 
of course be modified by changing the parameters (see Section 0) . 

Complementary results for SNR curves as a function of the site for different temperatures 
are shown in the second panel in Fig. |^. The resonance effects are evidenced by the crossings 
of the different curves. The crossings reflect the rise and subsequent drop in the SNR at 
a given site, and the fact that the curves cross at different sites confirms that the thermal 
resonance temperature varies from site to site. 

The non-monotonic behavior of the SNR with temperature is due to the fact that while 
the thermal power increases with T, the signal power at first also increases, and more rapidly 
than the thermal power. This is exactly the condition for the existence of stochastic reso- 
nance formulated by Dykman et al. p6t^p3| for single potentials, but that theory has not 
yet been generalized to extended systems. The fast increase of the signal power is consistent 



with trends found earlier for transmission of pulses along discrete arrays [13|. We found that 



a pulse in a hard anharmonic chain travels more rapidly with increasing temperature (except 
for the lowest temperatures, see below), while in a harmonic array a change in temperature 
has no effect on pulse speed. While a pulse broadens with increasing temperature in both 
arrays, the pulse in the anharmonic array remains relatively more compact. At sufficiently 
high temperatures the signal power at a given site becomes independent of T (as in the har- 
monic chain) because the signal response has reached its maximum value; a further increase 
in temperature only affects the signal further down the chain. The "crossover" temperature 
of the signal power from the increasing to the saturated behavior depends on the chain and 
signal parameters, a dependence explored in the next section. 

The fact that the SNR for the sites shown is very low (essentially zero) for the lowest 
temperatures at sites beyond the first few shows that for a purely anharmonic chain the 
signal essentially stops beyond the first few sites, whereas in the harmonic example the signal 
reaches all the sites shown even at the lowest temperatures shown. This is also consistent 
with the behavior shown earlier for the transmission of a pulse [0: the pulse velocity at 
sufficiently low temperatures is actually lower in the hard chain than in the harmonic case. 
Had we included a harmonic potential contribution equal to that of the harmonic chain, 
the SNR vs ksT curves would start at the same values as in Fig. ^ but for sufficiently 
distant sites from the first they would still be non-monotonic. We have omitted a harmonic 
contribution to present the thermal resonance effect in its purest form. 

The apparently monotonic behavior of the first few sites is due to the fact that for the 
temperatures shown the signal reaches these sites in any case. A resonance at these sites 
would be seen for different parameter values and/or at even lower temperatures. 

Our preceding descriptions point to another interesting measure of a thermal resonance, 
namely, the propagation length A, defined as the number of sites (i.e., distance along the 
chain) for which the SNR exceeds a certain threshold value. Fig. ^ presents the temperature 
dependence of A for harmonic and anharmonic arrays and an arbitrarily chosen SNR- 
threshold value of —1.2. Since the harmonic array does not exhibit thermal resonance, a 
monotonic decay of A with increasing temperature is observed. On the other hand, the 
hard chain shows a maximum for a moderate temperature. Again, the particular values of 
optimal temperature {ksT f» 3) and optimal distance (A ^ 17) can be modified by choosing 
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different parameters and/or different SNR thresholds, but the quahtative behavior persists 
as seen in Fig. 
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FIG. 6. Propagation length A as a function of the temperature for the anharmonic chain 
(standard case) and harmonic chain with k = 1 (other parameters as in the standard anharmonic 
case). The threshold value is SNR = —1.2. 

V. PARAMETER DEPENDENCES 



The parameters that can be varied in our model are the amplitude A and frequency ljq 
of the velocity of the first site, the coupling parameter k', and the damping coefficient 7. 




FIG. 7. First panel: SNR curves for different sites (from j = 4 to 13, top to bottom) along 
the anharmonic chain as a function of the temperature for the lower-amplitude case: A = 0.25, 
/c' = 5, 7 = 0.2, ujQ = 1. Second panel: SNR curves for the weaker-coupling case: A = 0.5, k' = 3, 
7 = 0.2, iOo = l. 
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The first panel of Fig. |^ shows the SNR curves as a function of temperature for a lower 
amplitude than in the standard case. Comparison with the first panel of Fig. ^ shows that 
the overall SNR is now (of course) lower and that the resonance temperature at each site has 
increased. The latter behavior indicates that for lower amplitudes the crossover temperature 
from an increasing signal to a saturated signal increases as the signal weakens. The second 
panel shows results for weaker coupling. Again the overall SNR is lower and the resonance 
temperature higher at each site. A weaker coupling thus has the effect of weakening the 
signal. 




FIG. 8. SNR curves for sites 4 to 13 (top to bottom) for the higher damping case: A = 0.5, 
A;' = 5, 7 = 0.5, ujq = 1. Second panel: SNR curves for the higher frequency case: A = 0.5, k' = 5, 
7 = 0.2, LOO = 2. 

The first panel of Fig. ^ shows the SNR curves as a function of temperature for higher 
damping. Comparison with the first panel of Fig. |^ shows that the overall SNR is again 
lower than in the standard case and the resonance temperature higher. The same trends 
are observed with a higher driving frequency as shown in the second panel of Fig. |. Again, 
each of these changes leads to an effectively weakened signal. 

Figures ^ and || illustrate clear trends. Decreasing A or k', or increasing 7 or uq, lead 
to the following consequences: a) the SNR at any given site and temperature decreases; b) 
the SNR resonance temperature increases at any given site; c) at a given temperature the 
resonance occurs at a site closer to the first. These trends are consistent with those found 
earlier for transmission of pulses along similar arrays [|T^ . We found that a pulse in a hard 
anharmonic chain travels more rapidly and relatively more compactly with increasing ampli- 
tude, decreasing damping, increasing coupling, or increasing driving frequency. By contrast, 
in a harmonic array changes in temperature or in the signal and damping parameters have 
no effect on pulse speed. 

VI. CONCLUSION 

We have demonstrated new thermal resonances in simple one-dimensional arrays of 
masses connected by hard springs. Our chains are not overdamped and hence include inertial 
effects. We have shown that the distance and amplitude of propagation of a signal imposed 
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at one end of the chain can each be optimized by tuning the temperature of the system. 
The resonance behavior reflects the temperature dependence of the distance traveled by the 
front of a signal and the dispersion of this front once steady state has been reached. At a 
fixed temperature, any parameter change that leads to an increase in velocity of propagation 
tends to increase the range of sites where the signal can be detected as well as the SNR 
at a given site. Increased velocity of propagation is also associated with a lowering of the 
resonance temperature for a given site. 

A complete understanding of the behavior of these simple anharmonic arrays in a thermal 
environment requires and relies on a number of other inquiries, some of which we have 
undertaken. One concerns the distribution of energy, and the persistence and mobility 
of energy fluctuations, in thermal equilibrium |ll2|Jl^ . We have found that soft anharmonic 
chains (that is, chains with soft interaction potentials) experience greater energy fluctuations 
than harmonic chains, which in turn experience greater fluctuations than hard chains. This 
is a direct consequence of the virial theorem. We have also established that fluctuations 
are mobile in harmonic and hard chains but not in the soft. Most interestingly, thermal 
fluctuations travel most rapidly and remain localized over considerably greater distances in 
the hard chain. These results in turn can lead to very different transition rate statistics and 
effective damping coefficients for a bistable impurity embedded in each of these arrays at a 
given temperature \T^\. 



Another interesting line of inquiry concerns the propagation of an energy pulse along 
such arrays. As mentioned earlier, we have found that the propagation of an energy pulse 
in a hard anharmonic array can be enhanced by immersing the array in a thermal bath, 
and that hard anharmonicity in the springs causes a tight and persistent packing of the 
energy 



Almost all of our results are numerical. Analytic results for the harmonic array provide 
some insights and point toward the possibility that an approach based on linearization 
methods [^] together with a linear response theory might lead to analytic insights for 
nonlinear arrays. We are currently exploring these possibilities. 

The thermal resonances and other interesting behavior that we have found for simple 
hard anharmonic chains are likely to be prototypical and therefore applicable to many other 
discrete systems with anharmonic interactions. While many other systems have been inves- 
tigated in which a careful balance of parameters at zero temperature or a manipulation of 
external noise lead to interesting localization, resonance, and synchronization phenomena, 
our systems are among the simplest generic systems in which the temperature can be used 
as the tuning parameter to achieve similar effects. 
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APPENDIX A: HARMONIC CHAIN 



In this appendix we consider a harmonic chain and calculate analytically some of the 
numerical results presented or invoked in the body of the paper. Although the behavior of 
linear chains is in general well understood, it is useful to present results in the particular 
context of this work. These are not readily available in the literature. 



1. Zero Temperature, With Signal 

Consider first an infinite linear chain (—00 < n < 00) at zero temperature, that is, in 
the absence of fluctuations (but with the dissipative contribution). The equations of motion 
then are 

Xn = k{Xn+l + Xn-1 - 2x„) - 'JXn, (Al) 

with the signal xq = Asm^uot) applied at site n = 0. We conjecture the quasi-stationary 
solution 

= _ Ae-A^l"! cos(cjot + \n\b) + c„ (A2) 

where the superscript s denotes the presence of the signal. The constants /i and b are to be 
determined. The additive constant c„ must be independent of n - it simply represents an 
overall translation of the chain (because it is not anchored) but otherwise does not contribute 
to the velocity analysis. The velocities associated with Eq. (|A2|) are 



<(t) = Ae-I"!'^ sm{uot + \n\b). (A3) 



This result is clearly consistent with the imposed signal at n = 0. The solution (|A3|) can be 
used to calculate the zero-temperature signal according to Eq. (|^): 

/+00 
e-'^''^ {Xn{t)Xn{t + T)) dr . (A4) 
-00 

The average (which in this deterministic case is only over time) eliminates rapidly oscillating 
contributions, leaving only a (5-function-type of contribution. The result is the power spectral 
density 

5:(u;o) = 4^e-2|«l^ (A5) 



To find the constant fi we substitute (^) into (|A1|) and set the coefficients of e*"^"* and 



those of e equal to zero (one resulting equation is simply the complex conjugate of the 
other). This immediately leads to the relation 

- col = k{e-^+'^ + e'^"*'' - 2) - z7cuo. (A6) 

We make the substitution u = e~^+*^ and note the symmetry in u and 1/u: 
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k(u + ^-2j +00^- i'yuo = 0. (A7) 
Next we multiply through by u and solve the resulting quadratic equation: 



i'jtoo) ± J {lOq — 2k — i'-f cooy — Ak'^ 



2k 



(A8) 



Note that u^u^ = 1, so if one solution is e"'^"'"*^, then the other is e^"*''. Which is which is 
not clear at this point. ^ must be positive for a physically acceptable solution, and this will 
be used below to sort out the choice. 



To extract /i we calculate the ratio which is either e^^ or e If 



-fi+ib 



then u^/u^ = e~'^^. If, on the other hand, u+ = e^~*^ then u^/u^ = e^^. Which it is will 
be seen at the end when we inspect the magnitude of the result. 

In calculating the ratio Uj^/u*_ one must exercise caution in taking complex conjugates 
because of the square root in Eq. (|A8|) ; it is not appropriate to simply change every i to a 
—i in such an expression if one does not know the signs of the terms inside the square root. 
It is helpful to write each complex number in terms of an amplitude and a phase: 



{ul -2k- i-fuo) = 2kae''^ 



(A9) 



so that 



2ka 



{uj',-2ky + Yoo'o 



1/2 



3 = tan ^ f „ 



Similarly, with 



we find 



(cuo' -2k- i^uj^f - = APu^e^'^ 



(AlO) 



(All) 



4fcV^ = ujo uj'^iuj'^ -Ak- 7^)2 + 472(cu^ - 2kf 



1/2 



2e = tan' 



.1 -2^{uji-2k) 



uJo{ui-4k-j^)' 



Then 



and the ratio u+/u*_ is thus real. Multiplying top and bottom of this ratio by u. 

u+ U+U- 1 



(A12) 

(A13) 
gives 

(A14) 



Further and explicitly 
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|m_|2 = [-ae''^ - ue'') x {-ae''^ - z/e"^^) = ^ z/^ + 2az/cos(/? - e). (A15) 
All that remains is the evaluation of the resulting expressions. In summary 

+ z/2 + 2aiy cos{p - e) = e'^^ (A16) 

or 

+ + 2au cos{P - e) = 6+^^. (A17) 

The choice is determined by whether the result is < 1 (in which case the first equality holds) 
or > 1 (in which case it is the second). The quantities a, u, j3, and e are thus completely 
defined in terms of the system parameters. 

We have tested our numerical simulations against the prediction ( |A5| ) for large ranges of 
parameter values and have found agreement to five significant figures. 

The solutions presented in this section do not obey A^-periodic boundary conditions. It 
is fairly straightforward using an image-like method to construct explicitly periodic solution 
by applying the signal at 0, ±A^, ±2N, . . . but the changes would be exponentially small in A^ 
and unimportant for sufficiently long chains. The quasi-stationary behavior assumed above 
for sites not too distant from n = sets in long before the signal at n = reaches sites ±A^. 



2. Finite Temperature, No Signal 

Now we solve the equations of motion 

Xn = k{Xn+l + - 2Xn) - ^X^ + fn{t) (Al8) 

without a signal. We are again interested in the stationary behavior. 
We define the Fourier transform and its inverse, 

N-l 1 N-1 

y,= Y. ^-e^'^^^", = ^ E ^/.e-'™"^. (A19) 

n=0 (j=0 

Transforming Eq. ( |A18D immediately leads to 

y, + Ak sin^ + = F,{t) (A20) 

where the inhomogeneous term is the transform of the noise: 

N-l 

F,{t) ^ fn{t)e^^''''. (A21) 

n=0 

The solution of this second order inhomogeneous differential equation is of standard form. 
The initial conditions (which we will take to be zero) are unimportant since we seek the 
long-time behavior: 

y^it) = hm ^ fp.ir) (e-('')(*-) - e-(^)(*-)) dr (A22) 

" t^co [r2{q) - ri{q}\ Jo ^ ^ 
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and consequently 
ijlit) = lim 

where 



- ri{q)] Jo 



Ak sin^ ( ^] . 



(A23) 



(A24) 



The superscript is used to stress the absence of a signal. 

With this result and the correlation function for the thermal fluctuations that follows 
immediately from Eq. (^), 

{F,{t,)F,,{h)) = 27 ksT Sit, - h)N5,,.g,, (A25) 

we obtain upon integration and Fourier inversion the velocity correlation function 

C°(r) = (x„(t)x„(t + r)) 



N 



Hiq) - rj{q)] 



(ri(g)e"^('?)*-r2(g)e"2(«)* 



(A26) 



We note in passing that this reduces to the standard correlation function for a single har- 
monic oscillator in a heat bath, C°(r) ksTe''^'^, when the coupling coefficient A; ^ 0. 
The associated power spectral density is 



dr C°(r) cos(cur) 



N 



[r?(g)+a;2][ri(g)+a;2]- 



(A27) 



Figure |^ shows the predicted spectrum (smooth curves) and the corresponding simulation 
results for two different temperatures. The agreement is typical of broad parameter ranges. 



0.15 
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FIG. 9. Spectrum S^{lo) for any site at two temperatures, fc^T = 0.05 (lower curves) and 
ksT = 0.1 (upper curves). The smooth curves are the predictions in Eq. ( A27| ) and the jagged 
curves are the numerical simulation results averaged over 500 runs. The irregularities can be 
eliminated at great CPU time cost by refining the time and/or frequency meshes and/or increasing 
the number of runs. Other parameters: /c = 3, 7 = 0.5. 



3. Finite Temperature and Signal: SNR 

Finally, we deal with the full equation of motion 

Xn = k{xn+i + Xn^i - 2x„) - 7X„ + /„(t) (A28) 

with the signal Xq = Asm{ujQt) applied at site n = 0. The fluctuation /o is set to zero since 
the velocity of that site is fixed by the signal. 

The Fourier transform is defined as before, but in transforming the equations of motion 
we must take into account that (|A28|) does not hold for n = 0. We thus multiply ( |A28| ) by 



^2mqn sum ovsT 71 but Only from n = 1 to — 1: 



N-l 



x„e2-^" = y,- xo. (A29) 



n=l 



Proceeding in this manner we obtain in place of Eq. ( |A2CI| ) 



y, + Ak sin^ (^^ j + jy^ = U,{t) (A30) 

where the inhomogeneous term now is 

Uq{t) = Fg{t) + k{2xo - xi - xn-i) + ixo + xq. (A31) 

The terms containing xq and its derivatives are known; xi and xjy^i (they are equal by 
symmetry) would have to be found "self-consistently" ex-post (see below). 

Equation (|A3(]| ) is again a linear inhomogeneous differential equation with constant co- 
efficients. Its solution has two additive contributions: 

y,{t) = y'rit) + yT'it) (A32) 

The portion y^°"^{t) is the solution of the homogeneous equation, that is, with Uq{t) = 0, 
and the contribution y^^^{t) is due to the inhomogeneity. 

Suppose that Fq{t) were equal to zero. The solution in this case, which we denote by 
y^it) is exactly the solution of Sec. |A 1| ; that is, the inverse transform of 

ylit) = °™(t) + (A33) 



at long times must be precisely Eq. ([A^). Now consider the consequences of again including 
Fq{t). The homogeneous part of the solution has to be exactly as before, i.e., 2/^°™ = y^'^°"^ 
because the homogeneous part of the equation has not changed, the constraint at Xq also has 
not changed, and initial conditions are in any case immaterial at long times. However, the 
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inhomogeneous part of the solution does change, and it does so in two ways. First and most 
important is the direct contribution of the thermal fluctuations to the solution. Second, and 
less important, is the change in Xi and (which is the same as that of Xi) caused by the 
addition of the noise, and the change that this in turn imposes on the other displacements. 
This contribution is more difficult to calculate and affects only sites near n = 0. Ignoring 
the latter we thus assume the solution 



(A34) 



where y^^(t) is given in Eq. (|A33|) and in Eq. ( |A22| ). 

The correlation function and spectrum are then simply the sums of the correlation func- 
tions and spectra for the system with the signal but no fluctuations (zero temperature) and 
those for the thermalized system in the absence of the signal. Thus, the spectrum is the 
sum of ( [A 51 ) (appropriately weighted by a delta function that places it at ujq) and ( |A27D . In 
particular, the SNR for the harmonic system then finally is 



SNRnM = log 



10 



/ 



42 



2nfi 



\ 



TV 



N-l 



(A35) 



which decreases monotonically with T as well as with n. This result is consistent with the 
linear response theory approach introduced and widely apphed by Dykman et al. |Q. Com- 
parisons of this analytic result with numerical simulations are shown in Fig. |1^ (cf. Fig. 
Equation (|A35 ) clearly captures the correct behavior; the small discrepancies between ana- 
lytic and numerical results (which always lie above the analytic curves) refiect the omitted 
terms discussed above. 




-3 I ' ' ' ' ' ' ' ' ' 1 -4 I ' ' ' ' ' ' ' ' ' 1 

1 2 3 4 5 5 6 7 8 9 10 11 12 13 14 15 

j 

FIG. 10. First panel: Typical SNR curves as a function of temperature for sites 7, 8, 10, 11, 
and 12 along a harmonic chain with k = 3, 7 = 0.5, loq = 1.0, and A = 0.5. The numerical 
simulation results (thin lines) in each case lie above the analytic curves (thick lines). Second panel: 
SNR as a function of site for different temperatures. 
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